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Abstract 

A two-parameter family of discrete-time exactly-solvable exclusion pro- 
cesses on a one-dimensional lattice is introduced, which contains the asym- 
metric simple exclusion process and the drop-push model as particular 
cases. The process is rewritten in terms of boundary conditions, and the 
conditional probabilities are calculated using the Bethe-ansatz. This is 
the discrete-time version of the continuous-time processes already inves- 
tigated in [1-3]. The drift- and diffusion-rates of the particles are also 
calculated for the two-particle sector. 



PACS numbers: 05.40.-a, 02.50.Ga 



1 Introduction 



The asymmetric exclusion process and the problems related to it, including for 
example bipolymcrization [4]. dynamical models of interface growth [5], traffic 
models [6], the noisy Burgers equation [7], and the study of shocks [8,9], have 
been extensively studied in recent years. The dynamical properties of this model 
have been studied in [9,10]. As the results obtained by approaches like mean 
field are not reliable in one dimension, it is useful to introduce solvable models 
and analytic methods to extract exact physical results. 

The totally asymmetric simple exclusion model on a one-dimensional lattice 
is one of the simplest examples from which exact results can be obtained. Such 
systems consist of a lattice in which every site is either empty or occupied by 
a single particle. Particles can hop to the right, if their right-neighboring site 
is empty. The steady-state of such systems have been extensively studied, for 
continuous-time as well as discrete-time evolutions. Among the methods used to 
study the steady-state properties of such systems is the matrix-product ansatz, 
[11-15]. Various methods have also been used to study the time-dependent state 
of such systems. In [16-18], generalizations of the matrix-product ansatz have 
been used to study asymmetric exclusion processes. In [19], an explicit form for 
the conditional probability of finding particles on specific sites in a system of 
asymmetric exclusion process was obtained in terms of a determinant. 

In [19], the coordinate Bethe-ansatz is used to solve the asymmetric simple 
exclusion process on a one-dimensional lattice. In [1], a similar technique was 
used to solve the drop-push model, and a generalized one-parameter model 
interpolating between the asymmetric simple exclusion model and the drop-push 
model. In [2], this family was further generalized to a family of processes with 
arbitrary left- and right- diffusion rates. All of these models were lattice models. 
The behaviour of the latter model on a continuum was investigated in [3] . The 
discrete-time version of the asymmetric exclusion process was discussed in [23] . 
In [20 22] , a similar Bethe-ansatz approach was used to study exclusion systems 
consisting of several kinds of particles. 

Here we consider discrete-time asymmetric exclusion processes in a one- 
dimensional lattic-c;. The scheme of the paper is the following. In section 2, a 
system is introduced which consists of a one-dimensional lattice in which each 
of the sites are either empty or occupied by a single particle. A discrete-time 
evolution is introduced and it is shown that the interaction between particles can 
be substituted by a suitable boundary condition. In section 3, the conditional 
probability of occupied sites is obtained. In section 4, the drift rates for the 
two particle sector are calculated. In section 5, the diffusion rate for the two 
particle sector is calculated. Section 6 is devoted to the concluding remarks. 

It is seen that for large times, the results of the continuous-time evolution 
are recovered, namely that the drift rates tend to the no-interaction drift rates, 
while the diffusion rate is generally larger than the diffusion rate of the non- 
interacting system. 
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2 A family of discrete-time exclusion processes 
on a one-dimensional lattice 

Consider a one-dimensional lattice, in which each site is either empty or occupied 
by one particle. The probability that the first particle is in xi, the second 
particle is in X2, etc. is denoted by 

P{xi,X2, ■■■), xi < X2 < ■■ ■ 

The process is that each particle can hop to the right, with the probability a, 
if the its right-hand side neighbor is empty: 

A9 0A, with the probability a (1) 

Consider the following evolution equation and boundary condition for the two- 
particle sector. 

Pixi,X2,t+l) ^{l-a)^P{xi,X2,t) 

+ a{l-a) [P{xi -l,X2,t)+ P{xi,X2 - 1, t)] 

+ P{xi~\,X2-l,t), xi<X2, (2) 

and 

P{x,x) ^ XP{x,x + 1) + fiP{x - l,x), X + ti = l. (3) 

Eq. (|2Jl describes a system with a diffusion process which occurs simultaneously 
for all particles. This is in contrast to a system for which at each step only one 
particle can hop to the right (if its right-hand site is empty). In the latter case, 
terms proportional to would be omitted from the above equation. 
Using 0, it is seen that 

P{x,x + l,t+l) =[(1 + Aa(l - a)] P{x,x + l,t) 

+ a{l-a) P(x -l,x+l,t) 
+ [a'^ +fia{l-a)]P{x~l,x,t). (4) 

So it is seen that ij^Jl and Q describe a system where particles can push: 

AA9 MA, with the probability (3 (5) 

where 

(3 — {1 — a) + = a ^ Xa (1 — a). (6) 

One can use ^ and Q, to obtain pushing rates in multi-particle sectors as 
well. This is especially simple in two cases: (A = 1,^ = 0) and (A = 0, ^ = 1). 
In the first case, one obtains 

n 

P{x, x + l,...,x + n,t+l) =(1 - a) ^ a™ 

X P(. . . , X + m ~ 2, X + m, . . . , X + n,t) 

+ a"+iF(x-l,...,a; + n-l,i). (7) 
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It is seen that the rate of particles all hopping to right is simply the rate of 
one particle hopping to right, to the power of the the number of particles. This 
shows that there is no pushing. This is the simple exclusion process. 
For the second case, one obtains 

P{x, a; + 1, . . . , a; + n, i + 1) =(1 - q:)"+^ P{x, ...,x + n,t) 

n 

+ a^(l-a)"-™ 

m=0 

xP[. ..,x + m — l,x + m + l,...,x + n,t). 

(8) 

This shows that there is a pushing process, the probability of which does not 
depend on the length of the block: 

A • • • ^ • • • A, with the probabihty a. (9) 

This is the drop-push model. 

3 The conditional probability 

The n-particle analogue of 10), can be written as 

F(x,i + l) = (C/P)(x,i), 

= [(1 - a + aTi)---(l - a + aT„)P](x,i), (10) 

where 

{TjP){xi, ...,Xn,t) := P{xi, ...,Xj-l,.. .,Xn,t). (11) 

For the evolution equation (|10|l . the Bethe-ansatz solution (the eigenvector of 
U) corresponding to the eigenvalue u is 

u*(x) = [{1 - a + aTi) ■ ■ ■ {1 - a + aTnMix.), (12) 

subject to the condition 

4'(. . . ,Xj = X, Xj+i — X, . . .) =\ 5'(. . . ,Xj = X, Xj^i = a; + 1, . . . ) 

+ /i *(..., = X - 1, Xj+i = X,. ..). (13) 

Using the Bethe-ansatz 

vI/k(x)=^A,e'^^-^(''\ (14) 

where a runs over n-permutations and 

Ai = 1, (15) 
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one arrives at 

n 

u=Y[{l-a + ae-'''i), (16) 

and 

Aaaj = S{k„(^j^,k„^j+i)) A^, (17) 

where aj changes j to j + 1 and j + lto j, and leaves the other numbers between 
1 and n intact, and 

S{kuk2) = , (18) 

This derivation is essentially the same as that used in [1-3]. 

Using these, the conditional probability (of finding the particles at x at the 
t, when they have been at y at the 0) is obtained as 

/Ant. 
^*k(x)e-''-''«*, (19) 

where the integration runs from to 27r, for each of the fcj's. Also, to treat 
the singularity arising from S in one is supposed to multiply A and ji in the 
denominator by e""^. It is seen that the right-hand side is equal (5x,y (for x and 
y in the physical region). So, the right-hand side satisfies the appropriate initial 
condition and evolution equation for the conditional probability, and hence is 
the (unique) solution to the conditional probability. 

4 The drift rates 

In the two-particle sector, the one-particle probabilities are defined as 

Pi{x,t) := P{x,X2,t) 

X2>X 

P2(a:,t) ^ P(a;i,x,t). (20) 



Xl<X 

Prom these, 

Pi {x, t + 1) =(1 - a) Pi {x, t) + aPi{x-l,t) 

+ A a (1 - a) [P{x, x + l,t)- P{x - 1, x, t)], 

P2{x, t + 1) =(1 - a) P2{x, t)+a P2{x - 1, t) 

+ Ha{l-a) [P{x - 2, a; - 1, t) - P{x - 1, x, t)\. (21) 

Defining 

(Xi)(t) ■.= Y,xPi{x,t), 

X 

{X2){t):=Y,xP2{x,t), (22) 
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(the expectation value of the position of the first and second particles) one has 

{Xi){t + 1) ={Xi){t) +a-\a{l-a) i), 



where 



{X2){t + 1) ={X2){t) +a + ^ia{l-a) P,(l, t), 



Pr{x,t) :=^P(t/,j/ + x,t). 



(23) 
(24) 



Writing H19|l for the two-particle sector, 

d2fc 



P{xi,X2,t;yi,y2,0) 



(27r) 



i (ki xi+k2 X2) 



1 - A '^^ _ g-i fci 



„i (fc2 xi+ki X2) 



-i (ki yi+k2 2/2) 



1 - Ae*'=i - fie-^f"^ 
X (1 -a + ae-*'=0*(l-a + ae-''=^)*, 



(25) 



one arrives at 



P,(x,t) 



dfc 
2^ 



-ifc (1/2-1/1) 



X (l-a + ae-"=)*(l-a + ae**)*, 



(26) 



where, using H24|) . the summation over y is done, which leads to a delta function 
J(/ci + fc2), using which one of the integrations is carried out. 

A steepest descent calculation shows that if t is large and a; is not large, 



Pr{x,t) 



\J'Ka (\ — a)t 



(27) 



So, for large i, 

(Xi)(0 =(Xi)(0) + at~A 

{X2W) ={X2){0) + at + fi 
One also has 



^ , a (1 — a) t ^ . , 
2a/— ^ '— + C + o{l) 



2A/^^ii^ + C + o(l) 



(28) 



((X2) - (Xi))(t + 1) =((X2) - {X,)){t) + a (1 - a) P,(l, t), 
{X){t + 1) =(X)(<) + a + a (1 - a) P,(l, i) 

(m (Xi) + A (X2))(t + 1) {Xi) + A {X2)m + a, 



where 



W ■.= \{{x,) + {X2)), 



(29) 
(30) 
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is the expectation value of the position of the particles. So for all times 



ifi (Xi) + A {X2)m = {fi (X,) + A (X2))(0) + at, 
and for large times, 



(31) 



{{X,} {X,m - (Xi))(0) + 2J^ii-^ + C + o(l), 



{X){t)^{X){0)+at 
+ - A) 



a (1 — a) t 



+ C + o(l) 



(32) 



where C is a constant. So the drift rates are large t are 



Vi := 



d(Xi} 
dt 

dt 



=Q! — A 



a(l- 


a) 


nt 


a(l- 


a) 





(33) 



(Xi) and {X2) are the expectation values of the positions of the first and second 
particles, respectively, and Vi and V2 are their corresponding velocities. As t is 
discrete, these velocities are defined only when the X^'s are smooth functions of 
t, which happens at large times. 

The above equations show that the drift velocities of both particles approach 
a for large times. The reason is that at large times the particles are far from 
each other and effectively do not interact with each other. But the next leading 
terms in velocities are negative for the first particle and positive for the second 
particle, which is expected from the hindering effect of the second particle on 
the first, and the pushing effect of the first particle on the second. One can see 
that the results obtained in [3] are recovered, provided one replaces a (1 — a)t 
with t. 



5 The diffusion rate 

Starting from (|^ . and defining 

X 

{xim ■.=j2x^ P2{x,t), 



(34) 
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one has 



{Xf){t + 1) ={Xf){t) + 2a {X,){t) + a 



Xa{l - a)^{2x + 1) P{x, x + 1, t), 



{Xl){t + 1) ={Xl){t) + 2a {X2){t) + a 



/i a 



(1 - a) ^(2x + 3) P{x, X + l,t). 



(35) 



Defining 



{X') :=-((X2> + (X|)), 



A' ■.={X^) {{X))\ 
(A^ is the variance of the position of the particles) one arrives at 
A2(t+1) =A2(<)+a(l-a) 



(36) 



+ a(l-a)(Ai-A) 



\ ^(2x + 1) P{x, x + l,t)- P,(l, t) {X){t) 



+ a (1 - a) - a - A)] Pr{l,t) 



a^ {\ - af {^i - \f 



(37) 



From H27|l . it is seen that the last two terms on the right-hand side vanish as 
t ^ oo. One also has 



l)Pix,x + l,t)=J2 [ 0(2a; + l)e'('=i+'=^)" 

X 

1 - Ae**^^ - ^e-''^'! 



A ki 



1 — A e* '^i — /i e * '^2 
X (1 ~a + ae-''''Y (1 - a + a e^**==)*, 



-i (fci yi+k-i y-i) 



)] 



1 - \e""' 



1 — A e* '^i — /i e ' 
X [l-a + ae-'^^y [l-a + ae"^^Y. 



(fei i;i+fc2 y2) 



(38) 



Here 5' is the derivative of 5 with respect to its argument. We are seeking those 
terms on the right-hand side of (|37() . which don't vanish as t ^ oo. On the 
right-hand side of H38|l . it is seen that the terms proportional to 5{ki + ^2) in 
the integrand, give rise to terms proportional to t^^/^ (for large t). (In fact 
these terms are proportional to Pj.(l, t).) Denoting the terms coming from that 
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part of the integrand which is proportional to S'{ki + ^2) by /, one has 



/ = 



1- Xe'''^ 



X e 



1 - Ae'*^! - 



(39) 



Differentiation of the fraction, S{ki,k2), needs some care. The derivative of 
this fraction is singular at ki — k2 — 0. To remove this singularity, one has to 
replace ki with fci + i e in the denominator. This prescription guaranties that 
in (|25|) . the integral of the second term on the right-hand side tends to zero as 
X2 tends to infinity. We are interested in the behavior of / for large times, and 
this is determined by the behavior of the integrand multiplier of S{ki + fe) for 
fci — ~k2, where k2 is small. One has 



dk2 \ 1 — A '^'1 — /i e 



ki— — k2— — k^0 



1-e 
/X — A 



i k 



k^O 



I e 



fc~0 



= (M-A)pf (^-j +t7r{^i-X)S{k). 

(40) 

Here pf means a pseudo-function (the Cauchy principal value in integration). 
Putting this in (|^ . and keeping only terms which don't vanish as i ^ cxo, one 
arrives at 

1= I —[-t7rip~X)S{k) + 2i-iat){l-a + ae'''Yil-a + ae-'^Y] + --- 

TT 



<M-A) + 



2at 



y^TTa{l — a)t 



So, 



^(2x + 1) P{x, x + l,t) = X) 



2a t 



Using (|?7jl and one arrives at 



y^Tra{l — a)t 



a t 



71" ^7ra (1 ~ a)t 
Putting (021) and lO in ^7^, one arrives at 

from which. 



lim — - — — a (1 — a) 



l + {li-Xf 













-I) 







(41) 

(42) 

(43) 
(44) 
(45) 
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This diffusion rate, is again defined only at large times, when one can treat 
as a function of continuous time. It is seen that it is in agreement with that 
obtained in [3], provided one replaces a (1 — a) i with t. 

6 concluding remarks 

The main result of the paper was to introduce a discrete-time discrete-space 
model, solvable through the Bethe-ansatz method. The model contains a free 
parameter (say A) that for certain values reproduces the simple exclusion model 
and the drop-push model. The conditional probability of finding particles at 
different sites was obtained as a function of time, from which in principle one can 
derive any correlation function. The conditional probabilities were calculated for 
the general multi-particle sector. The drift- and diffusion-rates, however, were 
explicitly calculated only for the two-particle sector, and it was shown that the 
results agreed with those of the continuous-time system at large-times. There 
remains the question of performing similar calculations for the multi-particle 
sector. For large times, one can put forward the following arguments. For large 
times, only the behavior of the integrand in H19I) around k = is important, 
and it is seen that for A = /i, there is no pole in the scattering matrix S. In 
fact, S becomes one for A = /i and k = 0, which shows that for A = /i and for 
large times, the conditional probability takes the form 

P(x,t;y,0) =^Po['T(x),i;y,0], (46) 

CT 

where Pq is the conditional probability of a system consisting of free particles 
hopping to the right, which corresponds to the one obtained with S* = in (f11?|l . 
This shows that for A = /z, at large times the system behaves collectively as a 
collection of free particles. Hence the drift- and diffusion-rates should be a and 
a(l — a) respectively, which agrees with the particular case of the two-particles. 

For the general case A 7^ /i, at large times the particles will generally be far 
from each other. The interaction terms coming from the scattering matrix, are 
in the from of products of two-particle scattering matrices. So it is plausible that 
for large times and for calculating up to 2-point functions, one neglects more- 
than two-particle interactions and interactions between non-adjacent particles. 
Then, from the n\ terms in the Bethe-ansatz solution for the n-particle sector, 
there remains only n terms. The drift rate at large times is expected to remain a 
again. For the diffusion rate, one can argue that it should be a{l — a) (the free- 
particle value) plus a function of (/x— A) which vanishes at /i = A. The additional 
term comes from the interaction of the neighboring particles. So it is expected 
to be proportional to (A — /i)^. This shows that the drift- and diffusion-rates 
obtained for the two-particle sector at large times, serve as qualitative results 
for the multi-particle sector as well. 
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